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Abstract. In Kaluza-Klein model with toroidal extra dimensions, we obtain the 
metric coefficients in a weak field approximation for delta-shaped matter sources. 
These metric coefficients are applied to calculate the formulas for frequency 
shift, perihelion shift, deflection of light and parameterized post-Newtonian 
(PPN) parameters. In the leading order of approximation, the formula for 
frequency shift coincides with well known general relativity expression. However, 
for perihelion shift, light deflection and PPN parameter 7 we obtain formulas 
Di:rg/[{D - 2)a(l - e^)], (D - l)rg/[{D - 2)p] and 1/{D - 2) respectively, where 
D is a total number of spatial dimensions. These expressions demonstrate good 
agreement with experimental data only in the case of ordinary three-dimensional 
(D = 3) space. This result does not depend on the size of the extra dimensions. 
Therefore, in considered multidimensional Kaluza-Klein models the point-like 
masses cannot produce gravitational fleld which corresponds to the classical 
gravitational tests. 
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1. Introduction 

The idea of the muhidimensionahty of our Universe demanded by the theories of 
unification of the fundamental interactions is one of the most breathtaking ideas 
of theoretical physics. It takes its origin from the pioneering papers by Th.Kaluza 
and O.Klein [Ij and now the most self-consistent modern theories of unification 
such as superstrings, supergravity and M-theory are constructed in spacetime with 
extra dimensions Different aspects of the idea of the multidimensionality are 
intensively used in numerous modern articles. Therefore, it is very important to 
suggest experiments which can reveal the extra dimensions. For example, one of the 
aims of Large Hadronic Collider consists in detecting of Kaluza-Klein particles which 
correspond to excitations of the internal spaces (see e.g. [3]). On the other hand, if 
we can show that the existence of the extra dimensions is contrary to observations, 
then these theories are prohibited. This important problem is extensively discussed 
in recent scientific literature (see e.g. [4]-|10j). 

It is well known that classical gravitational tests such as frequency shift, perihelion 
shift, deflection of light and time delay of radar echoes (the Shapiro time delay effect) 
are crucial tests of any gravitational theory. For example, there is the significant 
discrepancy for Mercury between the measurement value of the perihelion shift and 
its calculated value using Newton's formalism [TT]. It indicates that non-relativistic 
Newton's theory of gravity is not complete. This problem was resolved with the 
help of general relativity which is in good agreement with observations. Similar 
situation happened with deflection of light 'TT . The Shapiro time delay effect is 
used to get an upper limit for the parameterized post- Newtonian parameter 7 |13) . 
Obviously, multidimensional gravitational theories should also be in concordance with 
these experimental data. To check it, the corresponding estimates were carried out 
in a number of papers. For example, in [S], it was investigated the well known 
multidimensional black hole solution [M] and the authors obtained a negative result. 
However, this result was clear from the very beginning because the solution pH! does 
not have non-relativistic Newtonian limit in the case of extra dimensions. Definitely, 
in solar system such solutions lead to results which are far from the experimental 
data. The 5-D soliton metrics |TS]-[T7] were explored in In papers [5] and [5], 

it was found the range of parameters for which classical gravitational tests for these 
metrics satisfy the observational values. The black string (see e.g. jTS]) is a particular 
limiting case of such solutions with a trivial metric coefficient for the extra dimension. 
However, it can be easily shown that such solutions do not correspond to the point-like 
matter sources. 

In 5-D non-factorizable brane world model, classical gravitational tests were 
investigated in ^9j. Here, the model contains one free parameter associated with 
the bulk Weyl tensor. For appropriate values of this parameter, the perihelion shift in 
this model does not contradict observations. Certainly, this result is of interest and it 
is necessary to examine carefully this model to verify the naturalness of the conditions 
imposed. 

In our paper we consider classical gravitational tests in Kaluza-Klein models 
(factorizable geometry) with an arbitrary number of spatial dimensions D > 3. We 
suppose that in the absence of gravitating masses the metric is a fiat one. Gravitating 
point-like masses (moving or at rest) perturb this metric and we consider these 
perturbations in a weak field approximation. In this approximation, we obtain the 
asymptotic form of the metric coefficients. Then we admit that, first, the extra 
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dimensions are compact and have the topology of tori and, second, gravitational 
potential far away from gravitating masses tends to non-relativistic Newtonian limit. 
In the case of a gravitating mass at rest, the obtained metric coefficients are used 
to calculate frequency shift, perihelion shift, deflection of light and parameterized 
post-Newtonian (PPN) parameters. We demonstrate that for the frequency shift type 
experiment it is hardly possible to observe the difference between the usual four- 
dimensional general relativity and multidimensional Kaluza-Klein models. However, 
the situation is quite different for perihelion shift, deflection of light and PPN 
parameters. In these cases we get formulas which generalize the corresponding ones 
in general relativity. We show that formulas for perihelion shift, deflection of light 
and PPN parameter 7 depend on a total number of spatial dimensions and they are 
in good agreement with observations only in ordinary three-dimensional space. It is 
important to note that this result does not depend explicitly on the size of the extra 
dimensionf[|. So, we cannot avoid the problem with classical gravitational tests in a 
limit of arbitrary small (but non-zero!) size of the extra dimensions. It is worth noting 
that in paper [9] the authors arrived at the same conclusions in spite of they use the 
different approach. 

Therefore, our results show that in considered multidimensional Kaluza-Klein 
models the point-like gravitating masses cannot produce gravitational field which 
corresponds to the classical gravitational tests. 

The paper is structured as follows. In section 2 we get the asymptotic metric 
coefficients in the weak field limit for the delta-shaped matter gravitating source. 
These metric coefficients are applied to calculate the formulas of frequency shift, 
perihelion shift, deflection of light and PPN parameters in section 3. The main results 
are summarized in the concluding section 4. 

2. Weak gravitational field approximation 

To start with, we consider the general form of the multidimensional metric: 

ds'^ = gikdx'dx'' = 500 {dx°)^ + 2goadx°dx°' + gapdx°'dx^ , (2.1) 

where the Latin indices z, fc ~ 0,1,...,!? and the Greek indices a, /3 = 1, . . . , _D. D 
is the total number of spatial dimensions. We make the natural assumption that 
in the case of the absence of matter sources the spacetime is Minkowski spacetime: 
500 = f^oo = 1, ffoa = ma = 0, gap = rjaf} = Sap. At the same time, the extra 
dimensions may have the topology of tori. In the presence of matter, the metric is not 
a Minkowskian one and we will investigate it in the weak field limit. This means that 
the gravitational field is weak and the velocities of the test bodies are small compared 
to the speed of light c. In this case the metric is only slightly perturbed from its flat 
spacetime value: 

gik ~ rjik + hik , (2.2) 

where hik are corrections of the order In particular, hoo = 2ip/c'^. Later we will 

demonstrate that (p is the non-relativistic gravitational potential. The same conclusion 
with respect to ip can be easily obtained from the comparison of the non-relativistic 
action of a test mass moving in a gravitational field with its relativistic action. To 

X In the leading order of approximation, our formulas do not depend on sizes of the extra dimensions. 
All correction terms, where the sizes of the extra dimensions appear, are exponentially suppressed. 
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get the other correction terms up to the same order 1/c^, we should consider the 
multidimensional Einstein equation 

i?.. - (r.. - ^,9..T^ , (2.3) 

where So = 27r^/^/r(_D/2) is the total solid angle (surface area of the {D — 1)- 
dimensional sphere of unit radius), Gd is the gravitational constant in the (2? = D+1)- 
dimensional spacetime. We are going to investigate the weak field approximation 
where gravitational field is generated by N moving point masses. Therefore, the 
energy-momentum tensor is 

where mp is the rest mass and is the radius vector of the p-th particle respectively. 
All radius vectors r and are Z?-dimensional, e.g. r — {x^,x'^, . . . , x^) where are 
coordinates in metric (|2.1[) . The rest mass density is 

N 

p = ^mpS{r - Tp) . (2.5) 

p=i 

2.1. Xjc? correction terms 

Obviously, to hold in the right hand side of (|2.3p the terms up to the order 1/c^, the 
components of energy-momentum tensor ()2.4|) are approximated as 

Too«pc^ To„«0, r„^«0 ^ T = Tl^pc^. (2.6) 

Taking into account that /li^ are of the order of 1/c'^, the covariant components 
of the Riemann and Ricci tensors 

1 / 92g„„ d'^gu d'^gu d'^gkr. 



R ' 

* ™ 2 \dx^dx^ ' dx^dx"^ dx^dx"^ dx^dx' ^ 

+ 9np (^kl^ln - ^km^tl) ' ^kni = g'^ Riklm (2.7) 

up to the same order read correspondingly: 

* ™ 2 ydx'^dx^ dx'-dx"^ dx^dx"^ dx^dx^ ) ' 

km 2^ ydx'^dx'- dx^dx™ dx^dx^ dx'-dx'- 

_lf d^hl ^ d^hj d^h\ ,,d^hkm\ ^ 

2 \dx'^dx^ dx™dx^ dx'-dx™- dx^dx^ ) 
where h], = ri^"^hmk- With the help of the gauge conditions 

^ ^h^~\h\6^\^Q. (2.10) 



dx^ V'"' 2 
the formula (j2.9p can be written in the form 

1 .,d''hk 
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Taking into account that the derivatives with respect to = ct are much smaller 
than the derivatives with respect to x", we obtain from (j2.1ip : 

Roa ~ -A/ioq, ~ -A/Iq^ (2-13) 

where A = S"^d'^/dx"dx^ is the Z3-dimensional Laplace operator. It is worth noting 
that for the condition (|2.10p up to the order holds 

^{hi-lh[6^y0 + O{l/c^), g| = + O(l/c3). (2.14) 

Therefore, keeping in the left hand and right hand sides of (j2.3p terms up to the 
order we obtain the following equations: 

2SdGv . , „ 
Aftoo = n — P, A/lOa=0, 

= ■ ^2 ^ P^o,0 , (2.15) 

where = [2(13 — 2)/{D — l)]Gx>. Substitution of /iqo = 2lp/(? into the above 
equation for ft-oo demonstrates that (p satisfies the D-dimensional Poisson equation: 

A(^ = S'z3Gpp. (2.16) 

Therefore, is the non-relativistic gravitational potential. From (I2.15P we obtain 

It can be easily seen that in this approximation spacial coordinates of the metric (|2.ip 
are the isotropic ones, i.e. the spacial part of the metric is conformally related to the 
Euclidean one. It is worth noting also that the relation ha^/hoo = [1/(1? — 2)](5q^ can 
be also obtained from the corresponding equations in papers |141 120] . 

2.2. l/c' and correction terms 

Now, we want to keep in metric (|2.ip the terms up to the order 1/c^. Because the 
coordinate = ct contains c, this means that in goo ^^nd goa we should keep correction 
terms up to the order 1 /c* and 1 /c^ respectively and to leave gap without changes in 
the form g^p ~ rjap + hap with /iq^ from (|2.17p . 

First, we investigate the energy- momentum tensor (j2.4p which we split into three 
expressions: 

T™ = f: ^^c^[(_l)^g]-i/2^5(r - r,) , (2.18) 

^ rdt 
= Y: m,c[{~l)^g]-'/\f2^{r - r,) , (2.19) 

p=i 

^ rdt 

T-P = ^ „^^l(-lfg]-y\y^'^Sir - r,) , (2.20) 
p=i 



ft-Oa = 0, hap = — -hQodap = — —SaP ■ (2-17) 
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where = dx^ /dt . From (|2.20p we obtain up to order 1 (in units c) the covariant 
components 



(2.21) 



N 

Ta/s ~ ^ mpVpaVpisS{r - r^) . 
p=i 

Thus, taking into account the prefactor 1/c^ in the right hand side of (|2.3p . these 
components can contribute to gap terms of the order of 1 /c^ which is not of interest 
for us. For Tq^ we find from ({2^91) : 

N 

Toa « - X! mpCVpaS{r - r^) . (2.22) 

p=i 

Hence, these components can give in goa terms of the order of which is of interest 
for us. Finally, for Tqq we get from (|2.18p : 

Too = goigokT''' w {goofT°° 



2tp 



N 

X rUpC 
p=i 

N / 

X ^ rupc^ 1 
P=i \ 

N 



(-1)^(1 + ^) (-1 



2ip\ 



1 2ip 



c 

-I -1/2 



D-2 c2 



D 



1 -1/2 



(5(r-rp) 



1 + ^1(1 + 



N 



2c2 



2r/ ^ \^ f3D-4: 1 2\ 

2^mpc 5(r-rp) + _^mp ( — — — + -w^ I 5(r - rp) , 

p— 1 p—i ^ 



/SI) -4 1 



1 2ip 
D -2~^ 

3L» - 4 ^ 
15-2 ^ ^ 2^2 



<5(r-rp)«5]mpc2 l + :^^ + ^)5(r-rp) 



(2.23) 



where (/3p is potential of gravitational field in a point with radius vector r^. At the 
moment, we do not care about the fact that ipp contains the infinite contribution of 
the p-th particle. Thus, up to order 1 we get 



N 



p=i 

N 



D 



1 



£1-2 



"Pp - o^p S{r - Tp 



(2.24) 



With the help of ((2?23| and ((l^M)) we obtain up to the order 1/c^: 

N 



2SdGv 1 ^ 
-too - 1=; 75oo-t 



SdGz 



N 



D-V 

fiD-A 
V D-2 



SdGx 



D 

2{D - 2) 



^mp<5(r-rp) 



Similarly, from (|2.22l) and (I2.24p we get up to the order 

N 



2SdGv I I 



D - 1 SpGy 
D-2 c3 



^nipVpaSir - Tp) . 



(2.25) 



(2.26) 
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Now, we shall work out the left hand side of (|2.3p up to an appropriate orders of 
1/c. As we wrote above, we are looking for corrections of the order of 1/c* and 1/c^ 
to the metric components goo goa, respectively. To this end, it is convenient to 
present gik as follows: 

gik ~ Vik + hik + fik , (2.27) 

where /oo and /qq are of the order of and respectively. Then, the Riemann 
curvature tensor (|2.7I) reads 

1 / d^him . d^hki d'^hii d'^hk, 



R, 



ktn 



I / g'i J. Pt^ t. . f ., pii 

2 
1 



2 V dx^dx^ dx^dx"^ dx^dx^^ dx^dx^ 
I 



ydx^dx^ 


dx^dx™- 


dx'^dx"^ 


dx'-dx'- J 








d^hkm 


\dx^dx^ 


dx'dx™ 


dx^dx™ 


dx^dx'- 


^ d Jim 


d'fki 




d^fkm ^ 



1 

2 \dx^dx^ dx^dx"^ dx'^dx"^ 
1 np f dK 



—ri 

4 ' \ dx 



dhni 


dhki\ 1 


( dhpi _|_ 


dh 


dhijji 


dx'' 


dx"- J ' 


\dx^ 


dx^ 


dxP 


r)h 


dhkjYi 


\ f dhpi 


_l_ dhpi 


dhii 


dx'^ 


dx^ J 


1 \ dx^ 


dx^ 


dxP 



1 np ( (^^nk " • -inn ^ "■-Kill 1 j " • "pi "■-pi " • "U 1 /n nQ\ 



A' \dx 

From this formula we obtain the Ricci tensor 



^ 2^ dxWx' + 2 2^ dx^dx' + 2 

^ a ivu ( d^fT-im . d^hki d^ha d^hk„ 



fd^h, 




d'^hki 


^ ydx'^dx' 


dx^dx"^ 


dhnk ^ 


dhni 


dhki\ 


dx' 


dx'' 


dx" J 


dhnk 






dx'' 


dx^ 



2 ' ' ■'^ \dx''dx' dx'-dx-^ dx'^dx'^ dx^dx' 
1 il np f 9hn k dhrd dhki \ ( dh„i dhr,„r dhi, 



4 \ dx' dx'' dx"- J \ dx™- dx^ dx^ 



dhpi ^ 


dh 


dx™- 


dx^ 


f dhpi 


^ dhpi 


V dx' 


dx^ 



_l ii np ( dhnk dh„m _ dhkm \ f dhpi dhpi _ dha , , . 

4^^ \dx"' dx'' dx" J \dx' dx' dxP ' ' ^ ' ' 

where we introduced the notations: 

/ d'^h.m . d'^hki d^hu 



dx'-dx™- dx''dx"' 



and 



d'hl ^ d'hl ^ d'hj ^ d^hj d% 
dx''dx^ dx'-dx"' dx''dxP dx^^dx"^ dx''dx"^ 



P _ a I d'^firn , d'^fki d'^fii 

^ km — ^ 



dx''dx' dx'dx™ dx''dx^' 

d'f?n , d^^i , d^fH, , avf d^n 



dx^'dx'^ dx^dx"' dx''dxl^ dx^dx"' dx''dx"^ 
Taking into account that h^ — h'^ = 0, we get 



(2.31) 



Q2f^o d'^h" 

^ 9(# " 9(# ^^-^^^ 
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which are of the order and 1/c'^, respectively. For the components Fkm we obtain 

P ~ O •'0 _ P p ^ •'0 P_ (r) OA\ 

which are defined up to the orders f /c"* and 1/c'^, respectively. To get these expressions, 
we use the gauge condition 

dxP 2 dx^ ' 
Therefore, the 00-component of the Ricci tensor reads 



(2.35) 



2 2d{x^Y 2d{x^f 2d{x^Y 

1 1 92/1^ 1 .. , ^121. 



2 ' 2 9(a;0)2 ' 2 ' ' ^^^x^fe' 

^ V 52;' ax" y V dx^ dxP 

With the help of the following relations (which are correct up to the order 1/c^): 



92/100 


_ 1 


'"'dx^dx^ 


2 


f dhpQ 




V 





i?oo « + + • ^^A^ - ^(V(p)2 . (2.40) 



i'^ 9^j~2'' 9^9^^ ^^-^^^ 

and 

"■'""'1^ (^ - If - IJ) » If ^ - . , i.m 

where the condition (|2.14p was used in the latter expression, the OO-component (|2.36p 

of the Ricci tensor takes the form 

1 . 1 . , 1 2 . 2 

-A/00 + ^-^ ■ 

The Oa component of the Ricci tensor (|2.29l) up to the order 1 /c^ reads 

^ -2'' + 2^°" " 2^^ + 2^°" ^ 2^^"" + 2c^ 9^ ' 

where we used formulas p.33p and (I2.34p . 

Now, we come back to Einstein equation (12. 3p . Substituting (12.251) and (|2.40p 
into (|2.3p and taking into account (|2.5p and (|2.16p . we get the following equation for 
foo- 

A/00 + 4^^^ - 

2SDGvy^ f3D^4 D 
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With the help of the auxihary equation: 

4(V(^)^ = 2A((p2) - itpAtp (2.43) 
and equations (j2.5p and (|2.16p . equation (I2.42p takes the form: 

A (/oo -^y)^ ^ E ™p + 2(1^^2) ^0 - ■ ^'-''^ 

Here, tp'p is the potential of the gravitational field in a point with radius vector 
produced by all particles, except of the p-th. Substraction of the infinite contribution 
of the gravitational field of the p-th particle corresponds to a renormalization of its 
mass (see [21]). The solution of (|2.44p is: 

N N 

/oo = ^v'{r) + ^Y1 ^>'(r - r,) + • 1 " ^p) ' (2-45) 

p— 1 p— 1 

where (^'(r — r^) is the potential of the gravitational field of the p-th particle which 
satisfies the Poisson equation: 

= '^"'9^ = 5z,Gx,m/(r - r,) . (2.46) 

It can be easily verified with the help of (|2.5p and (|2.16p that (^'(r — r^) satisfies the 
condition: 

N 

<p(r) = E^'(r-rp). (2.47) 

p=i 

Therefore, substituting /loo — 2Lpj(? and /oo into (|2.27p . we obtain goo up to the 
order Xjc^: 

2 2 
5oo«l + — 'yj(r) + — (^^(r) 

p— 1 p— 1 

We should mention that the radius vectors Tp of the moving gravitating masses depend 
on time. In this case, potential (p{r) in (I2.47P also depends on time. 

The equation for /oq can be obtained by substitution of (j2.26p and p.4ip into 
Einstein equation p.3p : 

2{D ~ 1) SpGv ^ 19V ,o.Q^ 

p=i 

whose solution is: 

A.-^-ifi:w>-,)4ag.. (2.50) 

p=i 

where the function / satisfies equation 
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Therefore, substituting = and /qq into (|2.27p . we get ^qq up to the order 



2(D-1) 1 A , 1 d^f 

It is necessary to note that in the three-dimensional case D ~ 3 (|2.48l) and 
(|2.52p exactly coincide with (106.13) and (106.14) in 21 if we take into account that 
<f'{r - Tp) = -GNmp/\r - rp|. 

From now on we shall consider the case of one gravitating particle of mass mi = m 
at rest in our 3-D space but, for generality, moving with constant speed in extra 
dimensions. That is p = 1 ^ ip[ = and v" = dx°- /dt = (0, 0, 0, v^, v^^ vu), where 
V4^,v^, ...,vd are constants. In this case (|2.48p and (I2.52p are reduced correspondingly 
to 

300 « 1 + |^(r) + l^^(r) + ^^^f^^^{v) (2.53) 



and 



gap~-\l- 7^ — ^ • — V(i-) ) 5af) ■ (2.56) 



" ' (J - 2)c3 ^('•^ - ' (2.54) 

where 'y9(r) satisfies the Poisson equation: 

and = —gai3V°'v^ = w| -I- w| + ... + w|) + 0(l/c^) (at the same accuracy vp — —v^). 
Obviously, the transition to the case where the gravitating mass is at rest both in 
our three-dimensional space and in the extra dimensions corresponds to the limit 
Va — ^ = Q. In this case the potential ipir) as well as the function / do not 
depend on time t. We remind that the covariant components gap read (see (12.17^ ): 

1 2 
D~2 

To get all above results, wc did not use any concrete form of topology. The 
only things we used were assumptions of the flatness of metric in the absence of 
the gravitating masses and the weakness of the gravitational field and velocities of 
gravitating masses which perturb the flat metric. Now, to solve (j2.55p we should 
specify the topology of space and the boundary conditions. We suppose that the 
{D = 3 -|- (i)-dimensional space has the factorizable geometry of a product manifold 
Mo = R-^ X T'^. M.^ describes the three-dimensional flat external (our) space and 
T'^ is a torus which corresponds to a d-dimensional internal space with volume Vd- 
For this topology, and with the boundary condition that at infinitely large distances 
from the gravitating body the potential must go to the Newtonian expression, we 
can find the exact solution of the Poisson equation (|2.55p [23]. The boundary 
condition requires that the multidimensional G-p and Newtonian Gjq gravitational 
constants are connected by the following condition: SoGv/Vd — 4:TtGn- Assuming 
that we consider the gravitational field of a gravitating mass m at distances much 
greater than periods of tori, we can restrict ourselves to the zero Kaluza-Klein mode. 
For example, this approximation is very well satisfied for the planets of the solar 
system because the inverse-square law experiments show that the extra dimensions in 
Kaluza-Klein models should not exceed submillimeter scales [24] (see however [22] [23] 
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for models with smeared extra dimensions where Newton's law preserves its shape for 
arbitrary distances). Then, the gravitational potential reads 

/ A rgc^ 

where r^ is the length of a radius vector in three-dimensional space and we introduce 
three-dimensional Schwarzschild radius rg = IGf^mlf?. As we mentioned above, the 
gravitating mass m is at rest in our three-dimensional space but may move in the 
extra dimensions. In this case, the extra dimensional components of D-dimensional 
radius vector of the gravitating particle depend on time. The exact formulas for 
the non-relativistic gravitational potential (see [HI [53]) show that this dependence 
"nests" only in non-zero Kaluza-Klein modes which are exponentially suppressed in 
considered approximation. Therefore, in this approximation potential <p(r) in (j2.57p 
does not depend on time. 

It is worth noting that all the previous analysis works also in the case when the 
gravitating masses are uniformly smeared over some or all extra dimensions. Let us 
take for simplicity one (p = 1) gravitating mass m\ = m which is smeared over all 
extra dimensions. Obviously, this mass can move only in our usual three dimensions: 
v1 — dx^/dt = (w^, wj, uj, 0, . . . , 0) and its rest mass density (|2.5p now reads: 

p = l^m/ [] ,5(r3 - r(i)3) , (2.58) 

where Cq, are periods of tori. Then, the solution of the Poisson equation (12.16^ exactly 
coincides with the Newton potential if the multidimensional Gv and Newtonian Gm 
gravitational constants are connected as SdGv / Y^^^^i ~ ^t^Gn [22, 23 . Therefore, 
in this case the approximate formula (|2.57p becomes the exact equality: 

^(r) = ^(r3)^-^^-^. (2.59) 

2r3 

In the approximation (|2.57l) (or with (|2.59l) for "smeared" extra dimensions), the 
covariant components (|2.53|) . (|2.54p and (|2.56p take the form 

^ r3 2rl 2{D - 2)c^ r^ ' 

For the contravariant components we obtain: 

'I I rg 
r3 ' 2rl 2{D - 2)c^ r-i ' 

^ {D-2)c r^' ^ \ D-2 rsj ^ ^ ' 

It is not difficult to verify that these components satisfy the condition: 

_( l + 0(l/c^) + O(l/c5) \ 
9^k9 -(^o + 0(l/c5) 5o.p + 0{\/c^) ) ■ ^^-^^^ 

The metric components (I2.60p demonstrate that in this approximation the spacial 
section t = const is conformal to the Euclidean metric. Hence, the spacial coordinates 
are isotropic ones. It is convenient to use three-dimensional spherical coordinates 
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r3,0,ip instead of the Cartesian coordinates = x,x'^ = y,a;'^ = z. In these 
coordinates the metric (I2.ip reads: 



rg 2r2 2(L> - 2)c2 



(D-2)cr3 ^ 




Q=4 

1 r ■ 



i?-2 rg 



(dr^ + r^d6'2 sin^ 6ld7/>2 



- (l + ^) ((rf^^)^ + ('i-^)^ + . . • + {dx^f) ■ (2.63) 

As we mentioned above, this metric corresponds to a gravitating mass in the rest 
in our three dimensional space. If the mass is smeared over extra dimensions, the 
appropriate velocity components vanish. 

3. Classical gravitational tests 

Now, we want to check the obtained above multidimensional metric (|2.63p from the 
point of its consistency with the famous classical tests: frequency shift, perihelion 
shift, deflection of light and time delay of radar echoes (the Shapiro time delay effect). 
We also want to calculate the parameterized post-Newtonian (PPN) parameters for 
obtained metric coefficients. It is well known that four-dimensional general relativity 
is in good agreement with these experiments and observed PPN parameters. Can the 
considered Kaluza-Klein models with point-like sources also be in concordance with 
observations? 

3.1. Frequency shift 

To investigate the gravitational redshift formula in the spacetime ()2.63|) . we can use 
the famous expression for relation between the frequency ui of a light signal emitted 
at a point 1 with the metric component 500 1 1 and the frequency uj2 received at a point 
2 with the metric component gool2' 



(ffoo)^^^ =t^2 {gmY''' ■ (3.1) 



1 



a/2 



Therefore, up to the order 1/c we get 

c.2«^i(l + ^^) , (3.2) 

where non-relativistic potential Lp is given by (j2.57p . In considered approximation, 
this formula exactly coincides with the one from general relativity. Therefore, for this 
type of experiments it is hardly possible to observe the difference between the usual 
four-dimensional general relativity and multidimensional Kaluza-Klein models. 

3.2. Perihelion shift 

Let us consider now the motion of a test body of mass m' in the gravitational field 
described by metric (j2.63|) . The Hamilton- Jacobi equation 

g^^^^.rn'^^^Q (3.3) 
^ dx^ dx^ ^ ' 
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for this test body moving in the orbital plane 9 = n/2 reads 



l)w" rg dS dS 




(3.4) 



We investigate this equation by separation of variables considering the action in the 
form 



S ■ 



~E't + + SrAn) + ^"4 {x*) + - + Sd (a;-°) 



(3.5) 



Here, E' « m'c^ + i? is the energy of the test body, which includes the rest energy 
m'c^ and non-relativistic energy i?, and M is the angular momentum. Substituting 
this expression for the action S in the formula p.4p . we obtain an expression for 
{dSj-^/dr^Y holding there the terms up to the order 1/c^: 



/ dSr 



\ dra 

— r + 

_ US 4 



E' il- 



1 



1 



1 



D-2r3 
1 r„\^'^ 2{D - l)v' 



LA 



Dv^ 



2rj 



2{D - 2)c2 rs 



D-2r3 

( dSp 
\dxD 



dS 

{D - 2)c2 rg 'dx^ 
1 r„ 



1 



^-2 rg 



2m' E - {pI 



vl) 



E^ 

r,2 



Dm'^c rj 
2{D - 2) 



(3.6) 



,2 2 2{D-1) 
rn c r„ H z — m Ern 



1 

r3 



D 



D 



2{D - 2) 



,2 2 , 2(jJ-l) , ^ a \ 

' '■ff'' + {D-2) '^s 2^ ^ ) 



where Pa — dS/dx°' — dSa/dx" (a = 4, . . . , _D) are the components of momentum of 
the test body in the extra dimensions. If the gravitating and test masses are localized 
on the same brane then these components are equal to zero. Integrating the square 
root of this expression with respect to rg, we get in the following form: 

E^' 



2m' E 



{pI 



pI) 



1 



,2 2 2(1) - 1) , 
m c r„ -\ z — m Er,, 



D-2 



Dm''^c?r'jj 
2{D - 2) 



D ,2 2 , 2(1? -1) 

2{D^)" 



D > 

£1 = 4 / 



1/2 



drg 



(3.7) 



It is well known (see e.g. § 47 in ^3) that for any integral of motion / of a system 
with action S the following equation should hold: 



dS 

— — = const . 
ol 



(3.8) 
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Because the angular momentum AI is the integral of motion, the trajectory of the test 
body is defined by the equation 

+ const, (3.9) 



where we use (|3.5I) . 

Let now the Sun be the gravitating mass and the planets of the solar system be 
the test bodies. Then, the change of the angle during one revolution of a planet on 
an orbit is 

AV^ = -^A5.3, (3.10) 

where ASr^ is the corresponding change of 5*^3 . It is well known that the perihelion 
shift originates due to small relativistic correction e to Af^ in Sr^: M'^/r'l ^ 
(M^ - e)/rl. dSJl) shows that in our case e = Dm'^c^r'^/[2{D — 2)]. Expanding 
5*^3 in powers of this correction: 



(0) 



2M dM '■^ A{D-2)M 9M ' ^ ' ' 

where sif = 6*^3 (M^), we obtain 

Differentiating this equation with respect to M we get 

^-^"^-+ 2(i.-2)M2 ' ^'-''^ 

where we took into account —dASi^^^ /dM — A?/''^'^' — 27r. Therefore, the second term 
in p.l3p gives the required formula for the perihelion shift in our multidimensional 
case: 

Dnm'^c^rl DTTr„ 
'^=W^^^ iD-2)ail-e^) ^ ^'-''^ 

where in this equation we used the well-known relation = m'^rgC^a(l — e^) /2 with 
a and e being the semi-major axis and the eccentricity of the ellipse, respectively. For 
the three-dimensional case D = 3, this equation exactly coincides with formula (101.7) 
in [21]. It can be easily seen that the result (|3.14p does not depend on motion of the 
gravitating and test masses in the extra dimensions. 

It make sense to apply this formula to Mercury because in the solar system it 
has the most significant discrepancy between the measurement value of the perihelion 
shift and its calculated value using Newton's formalism. The observed discrepancy 
is 43.11 ± 0.21 arcsec per century. This missing value is usually explained by the 
relativistic effects of the form of p.l4p . However, only in three-dimensional case 
D = 3 (|3.14p gives the satisfactory result 42.94" which is within the measurement 
accuracy. For D — A and D = 9 models we obtain 28.63" and 18.40", respectively, 
which are very far from the observable value. 
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3.3. Deflection of light 



Let us consider now the propagation of light in gravitational field with metric (|2.63p . 
In the case of massless particles, the Hamilton- Jacobi equation p.3p is reduced to the 
eikonal equation: 



9 



= 0. 



which for the metric (|2.63l) reads 



(3.15) 




2{D - l)v°' rg (9* a* 
{D-2)c^ n 'dtd^ 



(3.16) 

where we take into account that light propagates in the orbital plane 6 = tt/2. The 
eikonal function ^I^ can be written in the form 

* = -LJot + + *,.3 (ra) + ^4 (x^) + ^5 (x^) + ... + *!) (x^) , (3.17) 

where ojq = —d'^/dt is the frequency of light and p is a constant. Later we will show 
that p is the impact parameter, i.e. distance of closest approach of the ray's path to 
the gravitating mass. Taking into account that k = ujq/c is the absolute value of the 
wave- vector, it is clear that M = pk = ploq/c plays the role of the angular momentum 
for the light beam. 

Now we consider the natural case when the light propagates in our three- 
dimensional space and does not have components of momentum in the extra 
dimensions, that is pa — d^'a/da;" =0, a = i, . . . , D. Then from p.l6p . using 
p.l7|) . we obtain up to the order 0(l/c^) the following formula: 

fd'fr.V 1 rA'W. r„ rl Dv^ 



V dr-s 

,2 



^ 1- 



^(1 



D-lrn 



D-2r3 



1 + ^ 
r3 



' g 
2rl 



2{D - 2)c2 rg 



? 9 

c^r^ 



D-2r3 



Integrating this expression we get: 



c 



1 



D-1 r„ 



1/2 



(3.18) 



(3.19) 



D-2r3 

Considering the term with rg/r^ as a small relativistic correction, we expand the 
integrand up to the order 0{l/c'^): 

D -I rgOjQ 



2{D~2\ 
D - 1 



f gWo ^ rg 



P ) drs 



2(D-2) c 
where the non-relativistic (i.e. gravity is absent: 

.2\ 1/2 

drs = 



arccosh — , 
P 



0) eikonal function is 



c 



1-^ 



1/2 



dr-i 



(3.20) 



(3.21) 
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For this non-relativistic approximation the trajectory of the hght beam is a straight 
hne. Indeed, in this case (by full analogy with p.9p ) we have 

r)\l>(") ,9\I>(°) 

W ^ ^^"^ + W = '^^"^ - ^rccosip/r,) = , (3.22) 

where the constant is taken in such a way that ^ 7r/2 for r3 — oo. Thus, the 
trajectory p = rscosi/j^"^ is the straight line. Obviously, in the non-relativistic case 
the total change of the angle is Ai/)^"^ = -dlS^f} jdM = tt. 

Coming back to the relativistic case p.20p . for the light beam travelling from 
some distance ra = i? to the closest approach to the gravitating mass at vs = p and 
again to the distance ra = R, the change of the eikonal function is 

A^rs ~ A*(") + !:^arccosh- . (3.23) 

D — 2 c p 

The corresponding change of the polar angle ij) is 

-KT7 =W+ -TTTT- = const 

dM ^ dM 



^ « + ^ !:£^ (i?^ - p^)-/^ . (3.24) 



Thus in the limit i? — ;> +oo we finally get: 

AV-«7r + §^^. (3.25) 
U — Z p 

Therefore, the second term in p.25p gives the required formula for the deflection of 
light in our multidimensional case: 

U — I p 

For the three-dimensional case D = 3, this equation exactly coincides with formula 
(101.9) in [21]. 

Now we apply this formula to the Sun. Obviously, the radius Rsun of the Sun is 
much greater than the size of the extra dimensions and approximation (j2.57l) works 
well on the distances r^ > Rsun- For general relativity and for a ray that grazes 
the Sun's limb (i.e. p k, Rsun) 5^ ~ 1.75 arcsec which is in very good agreement 
with observational data [H] . p.26p shows that we get this value of Sip only for usual 
three-dimensional space. In the case /) = 4 and D = 9 we obtain correspondingly 
Sijj w 1.31 and Sip ~ 1.00 which are very far from the observable value. 

3.4- Parameterized post-Newtonian Parameters and gravitational tests 

It is well known (see e.g. [251 HIl) that in PPN formalism the static, spherically 
symmetric metric in isotropic coordinates reads 

ds^ = f 1 - ^ + c'dt^ - (l + 7^) E {d. f . (3.27) 



1=1 



In general relativity, /3 = 7 = 1. However, simple comparison of equations (13.27^ and 
(|2.63p shows that the PPN parameters (3 and 7 in our case are 

/3 = 1 , 7 = ■ (3.28) 
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The latter expression shows that parameter 7 coincides with the corresponding value 
in general relativity ii D = 3. Only in this case 7=1. According to the experimental 
data, 7 should be very close to 1. The tightest constraint on 7 comes from the 
Shapiro time-delay experiment using the Cassini spacecraft: 7 — 1 — (2.1 ± 2.3) x 10~^ 
[l3[28l[29]. Ontheotherhand, for£) = 4,9wegetfrom(l3J8l)that7-l -1/2,-6/7 
respectively, which is very far from the experimental data. 

The formulas of the gravitational tests can be expressed via the PPN parameters 
[26j [28] . For example, the perihelion shift and the deflection of light read 
correspondingly: 

^^4(2 + 27-/3)^^, (3.29) 

d^j - (1 + 7) ^ . (3.30) 
P 

Now, if we substitute in these expressions the values from equation p.28p for f3 and 
7, then we exactly restore our formulas (|3.14p and p.26p . 

It makes sense to present the expression for the time delay of radar echoes (the 
Shapiro time delay effect) via the PPN parameters. This effect consists in time 
difference of propagation of electromagnetic signals between two points (or for a round 
trip) in the curved and flat spaces. Usually, a signal transmits from the Earth through 
a region near the Sun to another planet or satellite and then reflects back to the Earth. 
If the planet (or satellite) is on the far side of the Sun from the Earth (superior 
conjunction), then the formula for the time delay reads [261 128) : 

dt= {l + j)-^\nl I, (3.31) 

where r Earth and rpianet are the distances from the Sun to the Earth and to the planet, 
respectively. If we put into this formula the value of 7 from equation p.28|) . we get 

c , D 1 1 / '^'^ Earth '^planet \ /r» oi-»\ 

Obviously, this formula coincides with the general relativity only for D — ?i (see e.g. 
[37] )• For s-ll others values of D the time delay differs from the general relativity by 
the factor of 0(1). 



4. Conclusion 



In our paper we investigated classical gravitational tests (frequency shift, perihelion 
shift, deflection of light and time delay of radar echoes) for multidimensional models 
with compact internal spaces in the form of tori. We supposed that in the absence of 
gravitating masses the metric is a flat one. Gravitating point-like masses (moving 
or at rest) perturb this metric and we considered these perturbations in a weak 
field approximation. In this approximation, we obtained the asymptotic form of the 
metric coefficients. Until this point we did not require the compactness of the extra 
dimensions. This approach is valid for any number of spatial dimensions > 3 
and generalizes well-known calculations |21) in four-dimensional space-time. Then, 
we admitted that, first, the extra dimensions are compact and have the topology of 
tori and, second, gravitational potential far away from gravitating masses tends to 
non-relativistic Newtonian limit. It gave us a possibility to specify the non-relativistic 
gravitational potential for considered models. In turn, it enabled us to specify the 
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metric coefficients. In the case of a gravitating delta-sliaped body at rest, we used these 
metric coefficients to calculate frequency shift, perihelion shift, deflection of hght and 
parameterized post-Newtonian parameters /3 and 7. With the help of PPN parameter 
7, we also obtain the formula for the time delay of radar echoes. We demonstrated that 
for the frequency shift type experiment it is hardly possible to observe the difference 
between the usual four-dimensional general relativity and multidimensional Kaluza- 
Klein models. However, the situation is quite different for perihelion shift, deflection 
of light and the Shapiro time delay effect. In these three cases we obtained formulas 
which generalize the corresponding ones in general relativity. We showed that all of 
these formulas depend on a total number of spatial dimensions D and they are in 
good agreement with observations only in ordinary three-dimensional space _D = 3. 
This result does not depend explicitly on the size of the extra dimensions. Therefore, 
it is impossible to avoid the problem with classical gravitational tests in a limit of 
arbitrary small sizes of the extra dimensions. 

Therefore, our results show that in considered multidimensional Kaluza-Klein 
models the point-like gravitating masses cannot produce gravitational field which 
corresponds to the classical gravitational tests. Moreover, it is not difficult to show 
(see our forthcoming paper), that similar problem arises in the case of a compact 
static spherically symmetric perfect fluid with the following conditions for the energy- 
momentum tensor: Too ^ 2oq,T„^ , a, l3 — 1, . . . , D. To avoid this problem, it is 
necessary to break a symmetry between our three usual spatial dimensions and the 
extra dimensions. The branes are among the most natural candidates for solving this 
problem. Our results work in favor of the brane-world models. It is of interest also 
to check models with a non-linear action f{R). However, to prove viability of these 
models it is necessary to perform the similar investigations. 
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